Lithiation of InSb and Cu2Sb : A Theoretical Investigation 
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In this work, the mechanism of Li insertion/intercalation in the anode materials InSb and Cu2Sb 
is investigated by means of the first-principles total-energy calculations. The electron localization 
function for the lithiated products of InSb are presented. Based on these results the change in the 
bonding character on lithiation is discussed. Further, the isomer shift for InSb and Cu2Sb and thier 
various lithiated products is reported. The average insertion/intercalation voltage and the volume 
expansion for transitions from InSb to LiaSb and Cu2Sb to LiaSb are calculated and found to be 
in good agreement with the experimental values. These findings help to resolve the uncertainity 
regarding the lithiation mechanism in InSb. 



I. INTRODUCTION 



Intermetallic compounds present an attractive alter- 
native to graphite as anode materials in Li ion in- 
sertion batteries due in particular to the high capac- 
ity, an acceptable rate capability and operating poten- 
tials well above the potential of metallic lithium. The 
work in this direction gained tempo after the announce- 
ment that tin-based compoundsi*^ could operate effec- 
tively as alloying hosts for Li ions. In search of new 
anode materials one of the desirable qualities is that 
the structural and volumetric distortion of the host an- 
ode material on lithiation are small. This could be 
achieved if there existed a strong structural relation'^-^-^ 
between the host and its lithiated products. Some of 
the compounds which show such structural relations are 
CugSngSia, inSbiSi^iiaJii^, CuaSb", MnSb, MnzSbi^, 
SnSbi^, GaSbii etc. These compounds can be divided 
into three main structure types, namely NiAs (CueSns, 
MnSb), ZnS (InSb, GaSb) and CuaSb (MnaSb). In these 
materials, one of the components (Cu, In, Mn, Ga) is less 
active and is extruded on lithiation, while the other com- 
ponent provides the structural skeleton for further lithia- 
tion. In a previous work^ the lithiation of Ty'-CueSns was 
studied. In the present work, we investigate one repre- 
sentative compound from each of the other two structure 
types. These are InSb which belongs to the zinc blende 
structure and Cu2Sb which is the prototype for the third 
class of materials. 

Since it is desirable for an ideal anode material to re- 
versibly insert/intercalate Li ions while maintaining the 
structural stability, the knowledge of the exact amount 
of Li accommodated within the structure before the ex- 
trusion of the less active component is extremely impor- 
tant. In this regard, recently the electrochemistry of the 
Li insertion in InSb was studied experimentally, and con- 
tradictory results were published: 1) Vaughey et al.^ ini- 
tially proposed that 2 Li atoms could be incorporated 
within the InSb structure with very small volume ex- 
pansion (5.6%) before In extrusion, resulting in excellent 
capacity 2) Johnson et al.,^ subsequently reported that 
In extrusion occurred at an earlier stage according to a 
general reaction: {x + 2/)Li -I- InSb Li2.+j,Ini_j,Sb(0 < 



2;<2;0<?/<l)in which a Lij^Ini_ySb zinc-blende 
framework played an important role in the Li insertion/In 
extrusion process, but they could not determine the pre- 
cise amount of Li that could be inserted before the onset 
of In extrusion; 3) On the other hand, Hewitt et alm^ 
claimed, that x reached a maximum value of 0.27 before 
In atoms were extruded to immediately form LisSb in a 
two-phase reaction, and that Li2Sb could be generated on 
subsequent electrochemical cycling; 4) From in situ X-ray 
absorption measurements, Kropf et al^ and Tostmann 
et al}^ provided supporting evidence that Lia;_|_j,Ini_ySb 
electrode compositions were generated during the initial 
discharge of Li/InSb cells and that 40% of the extruded 
In was not incorporated back into the Sb lattice during 
the subsequent charging of the cells. In the light of the 
uncertainty of the reaction process, there is a need to clar- 
ify the exact lithiation path. To this extent we perform 
total-energy calculations for studying the mechanism of 
Li insertion/intercalation in InSb. A similar pathway is 
also provided for Cu2Sb. We also present the theoreti- 
cally determined isomer shift (IS) for these materials and 
various lithiated products. 

The change in the bonding character on lithiation is an- 
other important information that can be extracted from 
the first-principles calculations. Such analysis regarding 
battery materials, has been done in the past using the 
total^ or differenc»i^*i& charge density. But it was noted 
by Bader et alw^ and Becke et al^ that the informa- 
tion regarding the bonding in the molecular systems is 
not fully contained in the density (p) but rather in its 
Laplacian (V^p) and in the kinetic energy density. These 
quantities are used to calculate the electron-localization- 
function (ELF)fi^ which varies between and 1, with 1 
corresponding to perfect localization and 1/2 correspond- 
ing to electron gas like behaviour. Thus ELF in a way 
quantifies the bonding between various atoms in a solid 
and hence is ideal for studying the change in the bonding 
nature as new bonds are formed and (or) old bonds are 
weakened or broken by the insertion of Li atoms. We 
present the ELF for some of the lithiated products of 
InSb in the present work. 

The paper is arranged in the following manner. In Sec- 
tion II, we present the details of the calculations. Sec- 
tion IIIA deals with the energetics of lithiating InSb and 
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Cu2Sb. To facilitate the study of the change in the bond- 
ing character on the hthiation of these compounds the IS 
and the ELF are given in Section IIIB. Section IV pro- 
vides a summary of our work. 

II. METHODOLOGY 

Total energy calculations are performed using the full 
potential linearized augmented planewave (FPLAPW) 
method. In this method, the unit cell is partitioned into 
non-overlapping muffin-tin spheres around the atomic 
sites and an interstitial region. In these two types of 
regions different basis sets are used. The scalar relativis- 
tic Kohn-Sham equations are solved in a self-consistent 
scheme. For the exchange-correlation potential we use 
the local density approximation (LDA). All the calcula- 
tions are converged in terms of basis functions as well as 
in the size of the fc-point mesh representing the Brillouin 
zone. 

In order to compare the total energies, calculated from 
first principles, for any two structures one has to opti- 
mize the structural parameters. These parameters are 
the atomic positions and the unit cell volume, which 
are varied until the global energy minimum is found. 
Thereby the experimental crystalline data is used as a 
starting point. For all the energy differences reported in 
the present work these optimized crystal structures are 
used. For these calculations the existing WIEN2k codei^ 
is used. 

The details for the calculation of the IS of ^^^Sb ha-\ 
been presented before.^'' In the present work, the san 
value for the nuclear specific parameter (a — —0.846 : 
0, 002aQ mm/s) is used to convert the contact charge dei 
sities to the IS. The effect of the muffin-tin sphere size c 
the IS has also been studied previously.^" It was foun 
that the IS is not affected by the choice of the muffin-ti 
radius. A value of 2.0ao has been used in the presei 
work for all the atoms. 

The details for the ELF can be found in Ref. 
These calculations are performed using the newly deve 
oped EXCITING code. Calculation of ELF in a mixe 
basis like LAPW requires the kinetic energy density 1 
be made smooth across the muffin-tin sphere. Thus thes 
calculations are performed using super-LAFWj^i whic 
requires the matching of the wave function and its fir: 
and second derivatives at the muffin-tin sphere bounc 
ary making calculations computationally more demand- 
ing but affordable because the size of the unit cells is 
relatively small. 

III. RESULTS AND DISCUSSION 

A. Lithiation of InSb and Cu2Sb 

InSb exists in the zinc blende structure with In atoms 
at (0,0,0) and the Sb atoms occupying four tetrahedral 



sites: (0.25,0.25,0.25), (0.75,0.25,0.75), (0.75,0.75,0.25) 
and (0.25,0.75,0.75). This kind of structural ar- 
rangement leaves four tetrahedral sites vacant namely 
(0.75,0.25,0.25), (0.25,0.25,0.75), (0.25,0.75,0.25) and 
(0.75,0.75,0.75). Along with this there is another void 
in the structure at (0.5,0.5,0.5). On lithiation of InSb 
(or other zinc blende compounds), there are two distinct 
probable processes: 1) The incoming Li ions go to one 
of these voids. 2) A substitutional reaction takes place, 
i.e.. In (or the less active component) is extruded and is 
replaced by the incoming Li ions. In light of the existing 
controversy regarding these processes in InSb, we studied 
various possible lithiation paths. Using first-principles 
total-energy calculations we determine the formation en- 
ergies for various lithiated products of InSb using the 
following formula: 

yLi + LixInSb ^ Lix+yIni_ySb -I- yin (1) 

The formation energy for the compound Lix+yIni_ySb is 
given as 

AE = ELi,^,I„,_,Sb + jKn' - [ELi.InSb + jT^T]- (2) 

Here Ec is the total energy of the compound C and E™°* 
is the energy per atom for metallic A. The formation en- 
ergy defined this way reflects the relative stability of the 
compound Lix+yIni_ySb with respect to the reactants Li 
and LivInSbi^ 




Li concentration 



FIG. 1: The formation energies in eV per formula unit. 
The compounds formed from the insertion/intercalation of 
Li atom in the previous compounds are represented by circles 
and compounds formed by the substitutional reaction are rep- 
resented by triangles. 

These formation energies for InSb and its lithiated 
products are presented in Fig. 1. The compounds formed 
via the insertion/intercalation reaction of Li with the 
previous compound are shown as circles, while the com- 
pounds formed via the substitutional reaction are marked 
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InSb to LisSb is 0.83V, which is higher than the experi- 
mental value of the plateau in the voltage profile (0.6V— 
and 0.7Y^°). 

Another important thing to be noted from Fig. 1 is the 
difference in the formation energies between the two con- 
secutive compounds. After the formation of Li5/4ln3/4Sb 
it is energetically less demanding for Li atoms to come 
in and (or) for In atoms to extrude from the structure. 
Also in the case of the reverse reaction it is energeti- 
cally more demanding for In atoms to reincorporate in 
the structure beyond a concentration of 75%. This is in 
accordance with the experimental findings of Tostmann 
et al-r^ that only about 60% of the In atoms go back into 
the structure on charging and that Li atoms are extracted 
more easily and rapidly from this In deficient phase of the 
kind Lio.6+xIno.4Sb leading to a good cycling capacity. 



FIG. 2: (a) A part of Li-In-Sb phase diagram, (b) Close i; 
of the energetically favourable InSb-LiSSb part of the phai 
diagram 



by triangles. As can be seen the formation energies of ce 
tain compounds are negative, i.e. these compounds ai 
stable with respect to the reactants. A part of the Li-Ii 
Sb phase diagram is shown in Fig. 2. Fig. 2(a) consists ( 
the Gibb's trangle which consists of the energetically st; 
ble compounds in transition from InSb to LiaSb, whic 
are shown in detail in Fig 2(b). We find that insertic 
of 1/2 Li atoms per InSb formula unit leads to a volun 
expansion of 5.3%. In atoms are extruded beyond th 
concentration of Li, leading to very little further volun 
expansion. These results show best agreement with the 
experimental conclusions of Tostmann et ali^ The maxi- 
mum number of intercalated Li atoms before In extrusion 
as reported by Hewitt et alm^ is almost half the value ob- 
tained in the present work. According to the energetics, 
the compound Li2lnSb is unstable, but if formed leads to 
an enormous volume expansion of 15.4%. The total the- 
oretical volume expansion on going from InSb to LiaSb 
is 5.6% which is slightly more than the experimentally 
reported value of 4.4%. It can be noted from Fig. 1 that, 
like Hewitt ei aZ.jlS, we find the formation of Li2Sb during 
the first cycle to be energetically hindered. But we also 
note that the cubic phase of this compound is lower in 
energy than the hexagonal phase proposed earlier— 

The initial incoming Li atoms have the freedom to 
occupy either the site (0.75,0.25,0.25) or (0.5,0.5,0.5) 
in the InSb framework to form Lio.25lnSb. We find 
that the probability for Li insertion at the position 
(0.75,0.25,0.25) is higher than at (0.5,0.5,0.5). The dif- 
ference in energy between the two structures, Lio.25lnSb 
with Li at (0.75,0.25,0.25) and Lio.25lnSb with Li at 
(0.5,0.5,0.5), is 0.05 eV per formula unit. Our calculated 
average intercalation voltage (AIV) for a transition from 
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FIG. 3: The formation energies in eV per formula unit. 
The compounds formed from the insertion/intercalation of 
Li atom in the previous compounds are represented by cir- 
cles, and compounds formed by the substitutional reaction 
are represented by triangles. 

Cu2Sb exists in a tetragonal structure with space 
group P4/nmm and 2 formulae units per unit cell. The 
Cu atoms occupy the sites (0.75,0.25,0), (0.25,0.75,0.0), 
(0.25,0.25,0.27) and (0.75,0.75,0.73) while the Sb atoms 
are situated at (0.25,0.25,0.70) and (0.75,0.75,0.3). There 
exist voids in the structure at (0.5,0.5,0.5), (0.25,0.75,0.5) 
and (0.75,0.25,0.5). The two reaction mechanisms inves- 
tigated in the present are, 1) the insertion/intercalation 
of Li atoms, and 2) the substitution of Cu by Li atoms. 
The formation energies for these reaction paths are pre- 
sented in Fig. 3. Similar to the previous case, the Cu 
atoms are extruded from the structure after insertion of 
1/2 Li atoms per Cu2Sb formula unit. The first incoming 
Li atoms occupies the site (0.25,0.75,0.5) in the unit cell. 
Further lithiation leads to replacement of the Cu atom 
at site (0.75,0.25,0.0) by a Li atom to form LiCu3/2Sb. 
The optimized structural parameters for LiCu3/2Sb show 
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FIG. 5: Isomer shift in mm/s for the compounds on the en- 
FIG. 4: (a) A part of Li-Cu2-Sb phase diagram, (b) Close up ergetically favourable lithiation reaction paths for InSb. 
of the energetically favourable Cu2Sb-Li3Sb part of the phase 
diagram 

them easily identifiable by Mossbauer spectroscopy. 



a volume expansion of 18% per formula unit compared 1 
Cu2Sb. More incoming Li atoms leads to the formatic 
of the Li2CuSb and a free Cu atom. Li2CuSb exists i 
a cubic structure with space group F43m. In this un 
cell the Sb atoms are present at the origin, the face cei 
ters and four of the tetrahedral sites i.e. (0.75,0.75,0.75 
(0.25,0.25,0.75), (0.75,0.25,0.25) and (0.25,0.75,0.25) o. 
cupied by the Cu atoms. The Li atoms are located at tl 
other four tetrahedral sites and at (0.5,0.5,0.5). The thi 
oretically optimized structure of Li2CuSb shows a volun 
expansion of 26.4% per formula unit compared to the in 
tial compound Cu2Sb. This is in good agreement wit 
the experimental value of 25.2%.^ 

Energetically favourable Cu2Sb to LisSb part of tl 
phase diagram is shown in Fig. 4. The calculated vah 
of the AIV for going from Cu2Sb to LiaSb is 0.97 V, whic 
is slightly higher than the experimental value of 0.82 V- 
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B. Isomer shift and electron localization function 

An excellent method of studying the change in the lo- 
cal environment and symmetry is the Mossbauer IS. We 
have calculated the -'^^^Sb IS of various lithiated products 
of InSb and the results are presented in Fig. 5. The cal- 
culated value for InSb is in good agreement with the ex- 
perimental data of -8.6 mm/s.^^ Insertion reaction of the 
Li atoms with the host makes the IS steadily more nega- 
tive with a sudden increase when the In is extruded from 
the structure. The IS for the energetically hindered com- 
pounds LilnSb and Li2lnSb are -10.54 mm/s and -11.54 
mm/s respectively. These values are largely negative as 
compared to the other lithiated products of InSb making 



FIG. 6: Isomer shift in mm/s for the compounds on the en- 
ergetically favourable lithiation reaction paths for Cu2Sb. 

It can be seen from Fig. 6, the IS for Cu2Sb, like 
InSb, shows a decrease on lithiation via the insertion re- 
action and the continuous increase when the Cu extrusion 
starts. The IS has been shown to depend linearly on the 
partial number of s electrons (Ns)^*? and the variation of 
the IS with the Li concentration can be explained very 
well using this fact. As more Li atoms insert/intercalate 
in the host material, a part of the charge lost by them 
is gained by the Sb atoms. This causes an increase in 
Ns and hence a more negative IS. On the other hand, on 
lithiation via the substitution reaction there is no excess 
charge. This causes the decrease in and hence the IS 



5 



becomes less negative. Thus the change in the IS with the 
Li atom concentration is a good signature of the reaction 
mechanism and hence we suggest such experiments to be 
performed to shed more light on the reaction path and 
for the unambiguous clarification of whether compounds 
like LilnSb and Li2lnSb are formed or not. 



\\ — InSb 




(0,0,0) (1/4,1/4,1/4) (1/2,1/2,1/2) (3/4,3/4,3/4) (1,1,1) 



FIG. 7: Electron localization function for InSb, Lio.sSb and 
LigSb. 

The change in the bonding character on lithiation 
of the anode material can be studied using the charge 
density plots which have been used for CugSng^ and 
Cu2Sb^'*. In the present work we use ELF for such an 
investigation of InSb. Fig. 8 shows the ELF for InSb, 
Lio.sInSb and LisSb along the (111) direction. A cova- 
lent character of bonding between In atom at (0,0,0) and 
Sb atom at (1/4,1/4,1/4) can be seen as a strong local- 
ization of charge between the two atoms. Since only the 
valence electrons participate in the bond formation, the 
core charge has been ignored for all the atoms. The high 
value of ELF at the Li atom site is due the core states 
being taken into consideration as local orbitals. Lithi- 
ation of InSb to form Lio.sInSb leads to weakening of 
the bonding between In and Sb. This manifests itself 
by a lowering of the ELF from 0.9 to 0.79 between the 
two atoms. The bond between Li and Sb atoms is ionic 
in nature. This ionocity is retained all the way upto the 
formation of LiaSb. The presence of delocalized electrons 
between atoms is indicative of the metallic nature of the 
compound. 

IV. SUMMARY 

To summarize, the two prototype structures, InSb and 
Cu2Sb, for the use as anode materials in Li-ion batteries 



are investigated in the present work. From first-principles 
total-energy calculations we conclude the following: 

• On lithiation of InSb during discharge, the In 
atoms are extruded from the structure after inser- 
tion/intercalation of half a Li atom per InSb for- 
mula unit. The formation of Lii/2lnSb leads to a 
volume expansion of 5.4% 

• It is much easier for In and Li atoms to come 
out of and go into the structure after formation of 
Li5/4ln3/4Sb which supports the good cycling ca- 
pacity. 

• The formation of the compounds LilnSb, Li2lnSb 
and Li2Sb is energetically hindered. In particular, 
the formation of LilnSb or Li2lnSb would lead to 
enormous volume expansions of 14.4% and 15.4%, 
respectively, which is much more than any experi- 
mental observation. 

• Our calculated AIV for the transition from InSb 
to LiaSb is in agreement with the experimental 
value of the plateau in the measured voltage pro- 
file. These findings make our results closest to that 
of Tostmann et al^ 

• In the case of Cu2Sb, the Cu atoms are extruded 
after insertion of half a Li atom per Cu2Sb formula 
unit. 

• The trends in the IS with increasing Li concentra- 
tion is an excellent signature of the lithiation reac- 
tion mechanism. The formation of compounds like 
LilnSb and Li2lnSb can be unambiguously deter- 
mined using Mossbauer spectroscopy, since these 
compounds have large negative IS values compared 
to any other lithiated product of InSb. 

• The ELF shows that lithiation of InSb results in 
the weakening of bonds between In and Sb atoms 
which would lead, upon further lithiation, to an 
extrusion of In atoms. The bonding changes from 
covalent for InSb to ionic in nature for LiaSb. 
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